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Abstract 

o 

We describe the microstructure, shape and dynamics of growth of a droplet of martensite nucleating in a parent 
austenite during a solid-solid transformation, using a Landau theory written in terms of conventional affine, elastic 
deformations and non-affine degrees of freedom. Non-affineness, <j>, serves as a source of strain incompatibility and 
screens long-ranged elastic interactions. It is produced wherever the local stress exceeds a threshold and anneals 
diffusively thereafter. A description in terms of cj> is inevitable when the separation between defect pairs, possibly 
generated during the course of the transformation, is small. Using a variational calculation, we find three types 
. . of stable solutions (I, II and III) for the structure of the product droplet depending on the scaled mobilities of cf> 

parallel and perpendicular to the parent-product interface and the stress threshold. In I, is vanishingly small, II 
k> involves large <j) localized in regions of high stress within the parent-product interface and III where <f> completely 

wets the parent-product interface. While width I and size W of the twins follows I oc \fW in solution I, this 
relation does not hold for II or III. We obtain a dynamical phase diagram featuring these solutions and argue that 
they represent specific microstructures such as twinned or dislocated martensites. 
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1 Introduction 



Unlike the problem of solidification (Cahn and Haasen 19961, where a crystal nucleates and grows from a melt 



nucleation of a solid within another solid may encounter specific complications arising from crystallographic in- 
compatibility of the parent and product solids. In most cases it is geometrically impossible to fit a nucleus of 
the homogeneous product within the parent (austenite) matrix without gaps. To maintain continuity therefore, 
the system needs to adapt, which it does using several viable strategies. A common possibility is the simulta- 
neous nucleation of several crystallographically equivalent and degenerate variants of the product, producing a 



heterogeneous microstructure known as twinned martensite (Bhattacharya 20031, containing alternate arrays of 



elastically coupled variants or twins. At the martensite - austenite interface, crystallographic incompatibility is 
accommodated on the average. This is one of the various possibilities of accommodating incompatibility and there 



is an extensive literature on experiments (Cahn and Haasen 1996 



microscopic and coarse graine d theory jBarsch et al.||1984||Jacobs| 1985||Horovitz et al.| 1991| Rao and Sengupta] 



1997; Rasmussen et al. |2001||Rao and Sengupta 2003 Lookman et al. 2003 Paul et al. 2008| |Porta et al" 



Bhattacharya 2003 Olson and Owen 19921 



and coarse-grained and atomistic simulations ( Sengupta and Rao 2003 Rao and Sengupta| |2003| |Bhattacharya| 



[et~aL| [2008] |Kastner| [2003] [2006] |Sengupta et al.||2011[ ) of such martensites. 

On the other hand, there exists many other kinds of martensitic microstructures that incorporate defects, 



such as dislocation arrays (Olson and Cohen 19861, in order to make the interface compatible on the average. 



The incompatibility problem may also be resolved by large scale diffusion of particles, as is the case in ferritic 
microstructures which involve no twins and are produced at higher temperatures where solid-state diffusion is ap- 



preciable ( Cahn and Haasen 1996 Howell et al. 20111. In a realistic situation, multiple modes of accommodating 



interfacial strain compatibility may be present. Though there have been attempts to include the effect of isolated 
defects within an otherwise strain-only approach (Groger et al. 2008), a complete theory of microstructure, 



capable of describing dislocated, twinned, as well as ferritic microstructures has not yet been attempted. 

It is clear that one needs to incorporate additional degrees of freedom, other than strains, in order to describe 
multiple microstructures within a single framework. Accordingly, we include contributions from non-elastic, 
non-affine excitations (Kroner 1981| [Egami et aL 20071; i.e., local rearrangements of atoms which cannot be 



represented as a coarse-grained strain (Fig.|l|a)). Detailed molecular dynamics simulations of early time nucleation 
in a solid undergoing a (square — > rhombic) structural transition ( jBhattacharya et al. 2008[ l show that (1) non- 



affine strains are generated during the transition and are localized in regions we call non-affine zones (NAZ), (2) the 
transformation is driven by a small number of active particles (Sengupta et al. 2011 1 undergoing rearrangements 



within NAZs, (3) the NAZs lie close to regions where the non-order parameter (NOP) strains are larger than 
a threshold and lastly, (4) the dynamics of NAZ determine microstructure selection so that a martensite forms 
when the NAZs are localized at the interface while a ferritic microstructure forms when the NAZs spread over 
large scales. 

We had earlier described microstructure selection in solid state transformations, using a Landau theory which 
goes beyond strain-only approaches by including both elastic and non-affine degrees of freedom <j> that explicitly 
violate compatibility ( Paul et al.| 20081. We had shown that our theory could describe both martensitic and 
ferritic microstructures, and derive conditions when one or the other obtain. In this paper, we describe a droplet 



theory ( Rao and Sengupta 1997 Porta et al. 2009 I to study the structure and growth of a small region of 



martensite embedded in the matrix of the parent austenite phase. Simple but realistic variational ansatze for 
the order parameter and non-affine fields enables us to derive the optimized free energy of this droplet and show 
that there are multiple solutions corresponding to twinned and dislocated martensites within specific ranges of 
system-dependent parameters. 

The plan of the paper is as follows. In section[2] we discuss the dynamical equations for the elastic and non- 
affine variables. In section|3] we show that the equations of motion admit a Lyapunov functional, which we use 
to explore the space of solutions. In section[4]we describe the droplet variational ansatz of martensite embedded 
in the austenite phase. In section[5] we show that, in the space of our system-dependent external parameters 
the Lyapunov functional has multiple minima which may be related to structurally and dynamically distinct 
microstructures. We discuss our results in section[6]in the light of known results and conclude. 
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Figure 1: (a) Schematic illustration of a typical local rearrangement producing non-affine displacements. Note that 
such a transformation produces 5 and 7 coordinated neighbourhoods which may be interpreted as a pair of tightly 
packed dislocations with overlapping cores, (b) Illustration of the threshold dynamics for the non-afhne field given 
in Q. (c) Schematic illustration of NAZ produced when the stress &i{x) crosses the threshold o\ c . Horizontal lines 
show three values for er lc . For large a lc > cr\{x) (thin dashed line) everywhere, NAZ are never produced. As <t 1c 
decreases, small isolated NAZs begin to form — the regime described in our present calculation (bold dashed line). 
For (Ti c which is too small (dotted line), NAZs form everywhere and our droplet ansatz eventually breaks down. 



2 Dynamical equations for elastic and non-affine variables 

In this section, we will paraphrase our earlier derivation ( |Paul et al.| |2008[ ) of the dynamical equations for the 
elastic strains and (f> during a structural transition. While our theory is general, motivated by our earlier work 
phattacharya et al.[|2008| |Paul et"aL|[2008l|Sengupta et aT|[20TT| , we describe a formulation based on a square 
to rhombic transformation. The square to rhombic transition is a special case of the more general transition from 
square to oblique two dimensional Bravais lattices. The affine shear strain, e$ — du x /dy + du y /dx and the affine 
deviatoric strain e-2 — du x / dx — du y / dy are the order parameters (OP) for the square to oblique transition, where 
Ui(i — x,y) is the displacement field describing deformations from the (square) reference lattice. The equilibrium 
value of e2 for the rhombic lattice is zero. The volumetric affine strain ef = du x /dx + du y /dy is a non-order 
parameter (NOP) strain which relaxes much faster than the OP strains ( |Paul et al.[ |2008[ ) and will be taken to 
be slaved to the latter. 

We incorporate non-affine fields in our description by assuming that the total NOP strain may be decomposed 
into an affine e\ and a non-affine part 0, i.e., ei = e^ + 4>. This decomposition is analogous (though not equivalent 
!) to the decomposition of total strain into elastic and plastic parts, routinely made in the continuum theories of 



plasticity (Lubarda 2002 Acharya et al. 20081. In contrast to most plasticity theories, 4>, is however generated 
in the volumetric sector, an unusual requirement because deformation in solids is known to remain elastic for large 
compressional or dilatational stresses with plasticity being produced mainly under shear deformation, C3. This 
requirement can be readily understood, if we remind ourselves of the following : (1) the shear strain e^, which is 
the broken symmetry OP strain, is large within the martensite droplet and decays to zero in the parent phase, (2) 
gradients of e-$ generate large NOP (e^) strain at the interface in order to satisfy elastic compatibility and (3) the 
combined presence of and ei makes the interfacial region susceptible to local rearrangements of particles, this 
generates NAZs ([Bhattachar ya et al.| |2008) in regions of high NOP strain present at the moving parent-product 
front. The topological rearrangements, now parametrized by <f>, occur when the local conjugate stress 0"i({e3j) 
exceeds an appropriate local threshold (Ji c ({e3}) or yield criterion. Note that both the conjugate stress and the 
threshold which determine (j> are finally junctionals of e3. 



The Lagrangian is given by. 

r = f 2 



£ = / 7y [(«x + u 2 y )] dxdy - J 7 , (I) 
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where the first term is the kinetic energy with p the mass density, which may be set to unity without loss 
of generality. The free energy functional J- depends on the structural transition being considered and for the 
two-dimensional square to rhombic transition, one may set, 



r \ I dv.iv 



ai(ef + 4>) 2 + a 2 el + a 3 el — \^b 3 e\ 



+ \dA + ci(Ve? f + c 2 (Ve 2 ) 2 + c 3 (Ve 3 f] . (2) 

The parameters a\,a,2 and a 3 are elastic constants appropriate for the square (parent) phase and ci,C2 and C3 
are related to strain correlation lengths. For the range of parameters considered, T has three minima at e 3 — 
representing the square, parent or austenite and e 3 = ±eo = ±(63 + yjb\ — 4a 3 d 3 /2d 3 ) 1 ^ 2 , the two degenerate 
variants of the product rhombic phase. The connection with external control parameters such as temperature (T) 
is, as usual, through the temperature dependence of the coefficients. 

To obtain the equation of motion in the displacement fields, we need to solve the Euler Lagrange equation, 

d dc _ ac _ _dn 

dt dui dui diii 



The Rayleigh dissipation functional (Landau and Lifshitz 19861 TZ = \ J [71 (e^) 2 + 72^2 + 73e|l dxdy, where 
the coefficients 71 , 72 and 73 are the corresponding "viscosities" . 

We take to be a slaved variable which reaches steady state much faster than the OP strains. The equations 
of motion for the components of the strain can then be written as, 

V ex = q\ 3 — — - V 4>, (4 
oxoy 

e 3 = V 2 (^a 3 e 3 - b 3 el + d 3 el - c 3 V 2 e 3 ~ 20!^^- +y 3 e 3 ] 



+ 2 



dxdy 

■8&y{ a ^ + ^ ^ 
where qi 3 = 2(a2 — a 3 )/(ai + 0,2). The OP strain e2 has been eliminated using a modified St. Venant condition, 

V 2 ef - ( — - — \ e 2 - 2 d 2g3 = ~V 2 (j> (6) 
\ 9a; 2 dy 2 J dxdy 

where the strain incompatibility is completely determined by gradients of cj>. To obtain a closed set of equations, 
we need to prescribe a dynamics for the non-affine field cj>. To arrive at this equation, we note that <f> incorporates 
two dynamical processes. First, <j) is produced at a rate given by , whenever o\ locally exceeds a threshold <ti c . 
We denote the nonlocal response kernel as K(t,t'), where t,t' are times. This is the source of NAZs (Figfljb)), 
and goes to reduce the local stress o\. Second, as the transformation proceeds, <j> produced at regions of large 
<7i (Figjljc)) at the austenite-martensite interface, needs to anneal diffusively with coefficients c' x and c' y . The 
martensite droplet, being a broken symmetry phase automatically defines an interface or habit plane along which 
elastic compatibility is restored on the average. In general, atomic mobilities are different in the directions parallel 
and perpendicular to this interface suggesting c' x 7^ c' y . These processes suggest the following phenomenological 
equation of motion for (j>, 



c 



f ^H{\cj 1 \-a lc )K{t,t')dt' 

J -oo ft l 



■ 2 

x V dx 2 I v \ dy 2 



(7) 



In the above equation, the local volumetric stress cti = ST/Se^ — ai(ei + 0), H{x) is a unit step function and 
we have chosen a simple threshold criterion with a constant yield stress <ti c , independent of any state variable 
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(FigJTJb)) . The phenomenological parameters hi,c' x and c' y in may, in principle, be obtained from careful 
computer simulations of detailed atomistic models. 

3 The Lyapunov functional 

The equations Q and |5| together with |7| describe the dynamics of e, ; (i = 1,3) and 4> completely. We are 
interested in the steady state solutions of these equations. While the equations of motion of the strains is derived 
from the free energy functional the non-affme field follows a phenomenological dynamics |7| which does not 
arise from this free energy. In this section we show that, nevertheless, one may be able to define a Lyapunov 
functional whose minimizers do correspond to steady states of Q, (JsJ) and |7| thereby reducing our dynamical 
problem to one of optimization. 

Using K(t,t') — 5(t — t') to simplify the numerics, the phenomenological dynamics of the cj> field becomes, 

* - £-<w-~ >+<S) + 4(0), < 8 > 

where H (a) is the Heaviside function. Equation Q needs to be solved under the boundary conditions (f>(t = 0) = 
and <j> — > at large distances and reduces to two steady state situations depending on the magnitude of o~\. For 
| ci | > o~ic we have, 

m(e£ + 4>) _ , (d 2 4>\ , , (d 2 



hi \ dx 2 ) v \dy 

while for Icil < o\ c we obtain the anisotropic diffusion equation, 



^ + c ;(0)^o. 



x \dx 2 J 



Both these solutions may be combined as a single equation, 

c * (S) + Cy (0) + ai{4 + ^ )/(A ' ai) = 

where /(A, <ji c , <ti) belongs to a one-parameter family of smooth functions which represents the Heaviside function 
in the limit A — > 0. For our computational convenience, we choose 

*(£w£) + «tf + «{ i -»(-£)}-° (i2 » 

where we have absorbed the constant hi into the definition of the dynamical parameters c x and c y quantifying 
the ratios of the times it takes for (j> to form to that it takes to spread in the x and y directions respectively. 
One can show by direct plotting that the chosen representation for the step function, though strictly accurate 
for small o\ c , reproduces expected qualitative behaviour even if o\ c is larger. We emphasize that this choice is 
not a small ai c approximation, in principle we could have chosen a different smooth representation for the step 
function, accurate over a larger range of o\ c . 

We now seek a Lyapunov functional whose minimizers would automatically generate the steady state solutions 
for ei and (/>. It is straightforward to show that this corresponds to, 

1 f ( I J.2 , „ A ,1 , n ( Hi 



T L = -J |ai(^ +2ef0) + 2^)(ffi o + (Ti)e" 



+ H^J +Cy \d y/ 

+ ... (13) 



^dxdy 
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where the dots represent terms which are independent of <f>. Collecting terms we obtain, 



2a lc/ . 

(cric + CTije "i' +c x 

ai 



Ox 



+ Cy 



(14) 



One can check that Tl is, indeed, locally positive definite, and the system relaxes to the local minima of Tl\ 
the basic requirements of a Lyapunov functional (de Groot and Mazur 1984 Khalil 1996). Note however that 
the Lyapunov functional is not unique. Indeed, in a closed system the irreversible entropy production rate may 



itself serve as an appropriate Lyapunov functional ( Nieto-Villar et al. 2003 1 



We can now readily verify that a functional minimization of Tl obtains the required steady state solutions for 
strains and (f>. Finally, note that in (141, the first term within the integral is a purely elastic contribution whereas 
the second term is nonaffine. Elasticity favors the production of <j> (with appropriate sign) once the threshold 
criterion is satisfied, since the elastic energy is reduced through stress relaxation. On the other hand, large c x 
and/or c y inhibits the formation of (nonuniform) (f> at the NAZ where stress is large. The balance of these two 
tendencies optimize the magnitude and the sign of <j> and, consequently, the nature of the microstructure. 



4 The Droplet Ansatz 

Our variational ansatz for ez(x,y) is shown in Fig[2|a). We consider a rectangular droplet of the martensitic 
(twinned rhombic) phase occupying a region R3 of dimension L x W and consisting of a linear array of AT twins 
with alternating values of the order parameter es = ±eo represented by alternating strips of black and white. The 
droplet is oriented such that the habit planes or the parent product interfaces coincide with the lines x = ±W/2 
and the twin interfaces yu = —L/2 + kLN ', for k = . . .TV. We assume periodic boundary conditions in both 
directions. All interfaces are assumed to be flat which we later show to be a reasonable assumption for most cases. 
The region Rl represents the parent austenite (square) phase e$ = which is unaffected by the presence of the 
martensite droplet. Within the region Rl of size ~ L/Af, the order-parameter e$ still vanishes but the presence 
of the droplet produces stress fields a\ which determines the structure and dynamics of the interface. In Figj2|b) 
we have shown a typical interface of width n between twins with e% = +eo and e3 = — eo. The points are from a 
numerical solution of the dynamical equations Q (Paul et al. 2008), and the line is a linear interpolation which 



we use for our droplet ansatz here. The interface between the parent and the product where the order parameter 
e3 interpolates between ±eo and is similarly represented by a linear ansatz with interfacial width £. Within 
the region R2, the order parameter e$ vanishes and hence non-linearities in the free energy functional may be 
neglected. 

Since the volumetric strain is not a broken symmetry variable, it relaxes rapidly to its steady state value 
determined by the configuration of e$. We determine e± by solving the partial differential equation (jZJ) in the 
region R2. Equation Q is essentially a Poisson equation with "electric potential" being produced by a "charge 
density" d 2 e%/dxdy with the boundary condition = at distances large from the droplet, i.e., region Rl. In 
region R2, d 2 / dxdy = everywhere except at the triple junctions where twin interfaces meet the martensite- 
austenite interface (shown in Fig[2|a)&(c)). This analogy is expected to be valid in the presence of small or 
moderate amounts of <fi which reduces the effect of the anisotropic long range elastic kernel. This corresponds to 
a situation when the threshold stress is relatively large together with small diffusion times so that <j> accumulated 
at the interface is never too large (FigJTJc)). When, on the other hand, generous amounts of (f> are produced, this 
relation no longer holds. Therefore, in the limit of small interfacial width, e\ is the electric potential in 2D due 
to a configuration of point charges of alternate sign ±\/27rA centered at the triple junctions with, 



1 

: queo 



2tv U n 



l + l)y/e+rf- (15) 



For large Af, the potential (e^ ), at any point (x,y) within R2&iR3, due to this charge configuration in absence 
of 4> may be written as 

e?(x,y) = Xcos(qy) (e~^ +w ^ - e^*"^ 2 !) (16) 
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Figure 2: (a) Schematic diagram showing our ansatz for corresponding to the martensite droplet. Periodic 
boundary conditions are assumed in both the x and y direction. There are three regions Rl, R2 and R3. The region 
Rl contains the parent austenite (e^ = 0). All strain and associated stress fields are zero in this region. The region 
R3 is of size L x W, and contains AT twins (6 shown). The twins consist of regions where = ±eo, designated by 
white and black, which are degenerate in free energy. The gray regions are interfaces with widths r/ and £. The region 
R2 has no but contains a non zero non-order parameter strain ej 4 (and associated stress a\) which is produced due 
to the strain inhomogeneity created at the droplet. o\, which alternates in sign, reaches maximum absolute values 
at the blue and red colored circles marking regions of maximum strain inhomogeneity. (b) The ansatz e$(y) vs y 
(solid line) and the numerical solution of the dynamical equation (filled squares). The interfacial width 77 is shown, 
(c) A close up of the triple junction where the two twin variants (white and black) meet at a austenite-martensite 
interface. The red ellipse with widths fi x and f3 y represents a typical NAZ. 



where q = nAf/L and we have used a Fourier decomposition limiting ourselves to the most significant term. 

Initially, the stress o\ = aief is also largest at these triple junctions and is most likely to cross the threshold 
eric there generating <f> (which decreases <n). Accordingly, we supplant our ansatz for the OP field with one for <f> 
viz. 

4> = fa [ e -<*+^) 2 /^ _ e -(«-w/a) a /a£l J2 e -(v+^/^) 2 /^ 2 y cos (gj / ) (17) 

k = -Af 

which peaks at every triple junction with anisotropic widths fj x and /3 V in the x and y directions respectivelyand 
an amplitude which changes sign following that of e\ . The form of the ansatz (shown schematically in Fig(2fc)) 
interpolates between a <f) which may either be strongly localized (small /3 x ,/3 y ) or one which completely wets the 
austenite -martensite interface and is modulated with a period matching that of e\ and ez. 
Substituting the OP and slaved NOP strains and <j> in Tl we obtain the droplet free energy 

7-2 

T L = hLW + 2X pp L + X tw WAr + 2na 1 \ 2 - 7 + F^fo, (3 X , p y )Af (18) 
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Figure 3: (a) Plot of the scaled width j3 y of the 4* distribution for II and III (open and filled circles, axis on the left) 
and 4>q (open and filled squares, axis on the right) as a function of c y for c x at the triple point shown in FigQa) for 
(Tic = -001. (b) Typical plots of LjM as a function of VP 1 / 2 scaled by the interfacial width £ for a set of parameters 
over which each solution is stable. Bold line: solution I with eti = 0.05, A = 0.1, c x = 0.5 and c y = 4; dashed line: 
solution II with a\ = 1, A = 1., c x = 0.1 and c y = 1; dotted line: solution III with a\ = 1 and A = 1., c x = 1 and 
c y — 1. Note that while solution I reproduces the behavior of ferro-elastic martensites showing strong elastic coupling 
between the twins, in solution II and III the twins are elastically decoupled. 



where the constant Tb = aze\j2 — 6360/4 + d^e^/Q multiplies the first term representing the contribution from the 
bulk of the droplet. The second and the third terms with coefficients E pp = 63 dx (dez/dx) 2 = 2eo\/czX s , 
the surface tension due to the parent-product and T, tw = ^ dy {de^/dy) 2 = Aeoy/c^Ls for each twin interface 
with I s — azeo/3 -— 6363/5 + dsel/7, represent the energies of the parent-product and twin interfaces. The fourth 
term gives the contribution of the fringing NOP field ef. The function F^(<j)o, /3 X , /3 y ) in the fifth term is the 
contribution of the non-afnne field, </>, to the Lyapunov functional Tl which we had to obtain by numerical 
integration, checking our results against analytical expressions obtained in the pathological but instructive limit 
aic — > (see Appendix). F^{(j>o, /3 X , f} y ) contains terms which are linear and quadratic in </>o. Unimportant 
contributions from gradients of e\ have also been ignored for simplicity. 

Equation (18 1 is minimized with respect to /3 X , /3 V , J\f and <j)o as well as the interfacial widths £ and r\ for a 
fixed choice of all the parameters and size of the droplet L and W . 



5 Steady state solutions and the dynamical phase diagram 

Minimization of Tl with respect to the interfacial widths immediately yields 77 = £ = eo yf 63 ~/X s giving A = 
2gi3eo / ' \pK . Since we know that interfacial widths in solids are typically of the order of a lattice spacing 
( |Bhattacharya| |2003| ), this correspondence sets the microscopic length scale of our continuum model. Using 
these values, we minimize Tl with respect to the rest of the free parameters, viz., f3 x , f3 y , 0o and N for a range of 
values for c x , c y and o\ c . Our main results may be summarized as follows. 

1. We obtain three distinct solutions or dynamical phases, labelled 1,11 and III which correspond to unique local 
minima of Tl in particular regions of parameter space. 

2. The three dynamical phases show a wide variation in physical properties. While solution I is essentially 
affine, solutions II and III involve non zero values of (j>- 
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3. 4>, if present, is always strongly localized perpendicular to the habit plane, i.e. along x so that j3 x < /3 y <C W 
for all three solutions. Twin interfaces, therefore, do not ever contain any (j>. 

4. In solution II <p is localized at regions of high stress concentration at the habit plane, while in III, <j> completely 
"wets" the habit plane. 

5. Individual twins in I are elastically coupled while they are decoupled in II and III. 

6. While I is identified as twinned, ferro-elastic martensite, we propose that II and III are dislocated martensites. 

We now explain each of these results in detail. The magnitudes of 0o and /3 y are controlled by the dynamical 
constants c y and c x and o\ c . For fixed ai c , solution I, which exists for large c x and c y , has negligible non-affineness 
((po ~ 0) except at a small region of size ~ £ at the triple junctions. In this case, a suitable coarse graining of the 
interface taking £ — > may be used to derive a re- normalized strain-only theory without (f> (Paul et al. 20081. 



Therefore, this solution most resembles ferro-elastic martensite domains describable completely within strain-only 



approaches (Rasmussen et al. 2001 Lookman et al. 2003 Porta et al. 2009) 



and 



for both II and III as a 



To compare the nature of the two remaining solutions we have plotted /3 y 
function of c y for fixed c x over a typical range of parameters in Fig[3|a). The sign of (po in both solutions II and 
III is negative showing that the non-affine field effectively decreases the NOP stress (oc A). Solution II comprises 
localized NAZ of large \(f>o \ at the triple junctions while III has smaller \<j>o\ which completely wets the habit plane. 
As c y increases, |<^o| — > for solution III though before it vanishes altogether, a first-order dynamical transition 
to I intervenes. The width of the <f> distribution along the habit plane, /3 y is relatively insensitive over the entire 
range of c y over which these solutions are stable. 

Increasing M leads to smaller NOP fringing fields e± at the austenite -martensite interface and consequently 
lowers the contribution of the elastic energy due to elastic coupling between the twins. However, large M introduces 
more twin interfaces and hence increases the interfacial energy of droplet. For vanishing <j>, the interplay between 
these two factors alone decides the optimum number of twins A/". For example, we obtain 



L_ 

M 



aiA 2 ■ 



(19) 



This relation, between the size of the twin L/M and the width W, has been verified experimentally and has been 
obtained in earlier calculations ( Horovitz et al.||1991 Rao and Sengupta 1997 Porta et al. 20091. Substituting 
the minimized M in ( 18 1, one obtains Tl as a sum of bulk (oc LW) and surface (oc L + W) terms with an effective 



surface tension given by £ e / / = E pp 

Including 4> changes the relation between the twin size and the width depending on the resulting optimum 
value for cj>o since <j> changes both the numerator and the denominator of ( 19 1, viz. 



L 

M 



2ir(Y, tw W + A4q + 1 
aiA 2 + 2ai A0o 



(20) 



Where A, and B are positive definite functions of j3 x , /3 y and L/M. For viz. for solution /, we can 

negl ect the term quadratic in <f>o, the small linear contributions to both the numerator and the denominator in 
(201 then keeps the scaling relation unaltered. This is shown in Fi^3^b) where we have plotted L/M vs. W 1 ^ 2 
in scaled coordinates by changing W in our ansatz for a fixed L and for sets of parameters where each of the 
phases are stable. As \(j>o\ increases, the solution for L/M is not straightforward because the numerator is itself 
a (complicated) function of L/M and one needs to solve ( |20[ | self consistently and numerically. The result, also 
shown in Fig|3^b), reveals that (f> decreases stress, since non-affineness screens the elastic interactions between 
twins. For both solutions II and III, the twins appear to be completely decoupled elastically suggesting that 0o ~ W 
canceling the W dependence of the numerator making the twin width L/M almost independent of W . 

We obtain a dynamical phase diagrams showing regions of stability (lowest Tl) for solutions I, II, and III in 
Figs[4ja) and (b) in the scaled parameters c x /W£ and c y /c x for fixed o\ c . Remarkably, in these variables, the 
dependence on W, the size of the droplet, may be scaled away. For large c x and c y the system finds it increasingly 
difficult to generate non-affinity and there is a transition from structures II or III (large |<^o|) to I (small \ cf>o\)- Solids 
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Figure 4: (a) Dynamical phase diagram in the space of the scaled parameters c x /W^ (see text) and c y /c x showing 
relative stabilities of solutions 1,11 and III (see text) for three values of the twin width W/^ = 2500 (□), 3600 (o) 
and 4900 (A). Note that in the scaled coordinates, the phase boundaries (shown as red/gray lines) for different W 
collapse on a single set of three curves. The threshold stress a\ c = .001 (b) Dynamical phase diagram for three values 
of threshold o\ c — 0, .001 and .0025 for W/£ = 2500. Increasing <j\ c stabilizes solution I depressing the triple point 
and eventually removing II and III altogether. For o\ c > 0.005 the triple point obtains for the unphysical c y < 0. The 
rest of parameters common for both (a) and (b) are L/£ = 10 5 , with al = 1., a2 = 0.4, a3 = 1.25, b3 = —5, c3 = 10~ 5 
and d3 = 5.. 



which form only twinned martensitic microstructures are expected to have dynamical parameters corresponding 
to this region of our phase diagram. For lower c x and c y , on the other hand, <f> is produced and solutions II and 
III are stable. In this region, as c x is increased, there is a jump in the values of both the localization parameter 
2/3 y N ' I ' L and cj>o as the system switches from solution II with localized (j> with large amplitude to solution III with 
a relatively delocalized, broad (in y) (j> distribution albeit with a smaller amplitude. During this transition from 
II to III, f3 x also decreases, though this change is not as much as in j3 y . 

Finally, if a\ c is increased, predictably, non-affineness decreases stabilizing the "strain-only" I solution. In this 
case, the triple point in the phase diagram shifts to lower c y values and the phase boundary between the I and 
II phases shifts to the left (Fig|4^b)) and eventually tends to unphysical negative values for c y and c x . For large 
(Tic, again, one recovers the results of strain-only theories where twinned martensites with (j> — 0, i.e. solution I, 
is the only stable dynamical phase. 

While the interpretation of solution I as a strain-only, ferroelastic, twinned martensitic droplet is clear, the 
nature of solutions II and III still needs to be established. Our calculation, shows that <f> decreases the strong 
elastic coupling between the twins so that they are, in effect, individual grains of the product phase surrounded 
by an interface containing regions of large non-afhnity where the St. Venant's incompatibility \7 2 (j) 7^ 0. A 
non-affine region, defined in our work as a fluctuation in the local coordination number, may be interpreted as a 
compact droplet of dislocations with overlapping cores. Individual dislocations are easier to nucleate within such 



a region ( Bhattacharya et al. 2008 1 compared to the rest of the solid due to the lower barriers to dislocation 



nucleation within a NAZ. Normally the velocity of an austenite- martensite interface is large (Bhattacharya 



2003 1 and therefore nucleation of dislocations at the interface does not take place because nucleation times are 



too long. However, at high temperatures where the parent matrix is soft, and the interfacial velocity small, 
dislocation nucleation can and does take place with the strain mismatch across the parent product interface being 
accommodated by a lattice invariant strain composed of dislocations rather than twins. Such dislocated, "massive" 
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Figure 5: Comparison of the free energy as a function of the scaled average Burgers vector density b£/L obtained 
from solution II (dotted line) and direct evaluation using an array of dislocations (bold line) for a droplet containing 
a single twin with c x and c y near the triple point values for <j\ c — .001. Note that the two curves yield comparable 
values for b. 



martensites are commonly seen, for example, in low carbon, soft, iron alloys ( Olson and Owen 1992 Olson and 



Cohen 19861. The movement of the product interface, or habit plane, is then accomplished by "glide" motion of 
the dislocations while movement of dislocations along the habit plane is due to the much slower "climb" motion. 
In order to investigate whether the presence of non-affine fields at the parent- product interface can be indeed 



and Sethna 



and plastic, , parts. Accordingly we identify, the incompatibility V x e 1 ' 
tensor which is given in component form, 

de p 

f _ Jijk mj 

axi 

where e ljfc is the unit antisymmetric tensor, the summation convention is implied and the indices run from 1 — 3; 
for plane strain conditions relevant here, we need to look only at terms corresponding to k — 3, viz. 



2008) 



related to interfaci al dislocations, we use the correspondence suggested by ( Acharya et al 

^— ± ^ V 

2006 1 where the total strain is decomposed as iT = e + t? i.e. consisting of a sum of afiine 



and JLimkumnerd 



with the dislocation density 



C32 



i 

<l.r 



representing the x and y components of the Burgers vector density. To obtain the total dislocation density, we 
need to coarse grain £31 and £32 over a region L x £. The y component vanishes by symmetry yielding finally the 
mean Burgers vector content b per unit length of the interface. From the results of the last section, we conclude 
that b w for solution I; the austenite- martensite interface in this case being virtually free of any dislocation 
content. This is not true for solutions II and III where we expect to obtain an optimum b related to the optimum 
profile for the <f> field i.e. cj>o, (3 X and j3 y . Since all of these quantities are almost independent of c x and c y within 
the range over which each solution is stable, it, therefore, makes sense to compare b with that obtained from 
an array of evenly spaced dislocations at the austenite- martensite interface each with Burgers vector b oriented 

Considering that creating each dislocation needs an expenditure 



1986) 



along the x direction ( Olson and Cohen 
of core energy E core = /4ir one can obtain the Burgers vector density by optimizing the Landau free energy, 



T d = 



die? + a 2 e 2 2 + a 3 e3 + ^b 3 e3 + -d 3 el + c 3 (Ve 3 ) 2 



dxdy + N d E cor 



(21) 



11 




0?001 0.01 0.1 1 
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Figure 6: (a) Schematic diagram showing the modulation of the habit plane, (b) Plot of the curvature of the 
Lyapunov free energy Ti, as a function of the scaled modulation wavenumber for solutions II (filled corcles) and III 
(open squares). A negative value of this quantity signifies instability. Note that the interface continues to remain 
planar with the curvature decreasing roughly as ~ l/9m- 



where the strains have contributions from the OP (e3(x,y)) and NOP (e^(x,y)) strains from the martensite 
droplet as well as from the array of interfacial dislocations with each dislocation (chosen to be at the origin) 
contributes, 

~di = ~ Dv { x i + y*y (22) 

ay (x 2 + y 2 ) 2 

= ^= Dx x2 ~y 2 (24) 
dy dx (x 2 + y 2 ) 2 

with D — a,3b/4n(l — v) and v — (3ai — 2a^) /2(3a\ + 0,3) the Poisson's ratio. 

A comparison of the free energy as a function of the dislocation density for these two calculations with a specific 
values for c x and c y for a droplet containing a single grain (iV = 1) of the II phase is shown in Fig|5] Taking the 
magnitude, b, of the Burgers vector as an adjustable parameter we see that the optimum for the two alternatives 
coincide if b = 1.3£, an entirely reasonable value. Changing c x and c y causes only minor changes in the value 
of b. We therefore identify this solutions as representing dislocated or massive martensite. The correspondence 
with dislocations also lends meaning to the coefficients c x and c y as dislocation mobilities in the glide and climb 
directions (perpendicular and parallel to the habit plane) respectively, although exact numerical correspondence 
with measured mobilities for any particular real system is beyond the scope of our theory. Immediately, we note 
that portions of the dynamical phase diagrams shown in Fig|4] showing a non-affine phase for c y /c x 3> 1 may 
be physically unrealizable for most materials since climb is usually much slower than glide except for very high 



temperatures where an un-twinned ferrite microstructure with non-planar interfaces is probably more stable ( Paul 



et al. 2008 \ 



Our theory, incorporating non-affine fields therefore offers a common framework within which both twinned 
and dislocated martensites in a variety of systems with different dynamical characteristics, each corresponding to 
physically realizable portions of our dynamical phase diagram Fig|4j may be described and discussed. 
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6 Discussion and conclusion 



In this paper we have introduced a coarse grained continuum model which extends strain-only theories of marten- 
sites by including non-elastic degrees of freedom capable of describing localized regions called NAZ within which 
particles are known to have non-affine displacements. In the final section of this paper we discuss our present 
results, taking up a few issues which examine our assumptions, illustrate the relevance of our findings and point 
out avenues for further work. 



6.1 Stability of the parent-product interface 



Our ansatz, for the shape of the martensite droplet is rectangular implying flat twin and parent- product inter- 
faces. How good is this assumption? While twin interfaces are usually flat because they minimize elastic energy, 
moving parent-product interfaces during growth of the product phase may undergo dynamical (Mullins-Sekerka) 
instabilities producing a variety of interesting structures and shapes, eg. dendrites (Mullins and Sekerka 19641. 



On the other hand, austenite- martensite interfaces in ferro-elastics remain planar despite the often large velocities 
involved. This is easy to understand once we realize that the Mullins-Sekerka instability is strongly suppressed at 
wavenumbers below the inverse capillary length oc E e // oc y/W in the limit W — > oo. Since solution I corresponds 
to this case, our ansatz remains trivially valid. In order to investigate the validity of our ansatz for solutions II 
and III, we calculate Tl numerically after imposing a modulation ~ A m sm(q m y) (Fig|6^a)) on the habit plane. 
In Figj6^ta) we show that the curvature of the free energy d 2 Thfd!\ 2 m ~ q„ 2 > proving that our ansatz is 
stable against such modulations. Nevertheless, dynamical instabilities of the Mullins-Sekerka type can still arise 
since <fi screens elastic interactions leading to a less than \fW dependence of the effective surface tension. Indeed 



austenite- ferrite interfaces which contain many defects (Howell et al. 2011 1 are hardly planar. This has also been 



observed in our own MD simulations (Bhattacharya et al. 20081 and elasto-plastic theory (Paul et al. 2008) 



6.2 Stick-slip motion 



Non-afnneness at the parent- product interface may produce dynamical phenomena which are qualitatively new. 
For small <f>, one obtains steady growth of the interface, with <j> being advected along with the growing interface. On 
the other hand, our threshold dynamics for yield, may give rise to intermittent dissipation into non-affine and affine 
sectors producing stick-slip dynamics of the austenite-martensite interface which we would like to investigate in the 
future. Such stick-slip motion has been known to be prevalent in many driven dissipative threshold systems, e.g., 



frictional sliding (Heslot et al. 19941, peeling of an adhesive tape (De et al. 20041, earthquake faults (Burridge 



and Knopoff| [l967|), Barkhausen noise i n ferromagnetic transitions ( jPerkovic et al.| |1995|), serrations in stress- 



strain curve (Portevin and Le. Chatelier 1923 1 and rate independent hysteresis ( Abeyaratne and Knowles 19971 



The system usually spends a large time in the stuck state during which stresses are build up and subsequent rapid 
release of stress facilitates the slip of the system. Indeed, even structural transitions in solids, including martensitic 
transformations (jPinsook and Acklandj 2000 Ngan and Truskinovsky 1999), are known to be accompanied by 



stick-slip motion ( Buldum and Ciraci 1997 1 of the interface. The most predominant view is that stick slip in 
these systems arises from the presence of compositional disorder (Bhattacharya 19991. However, experiments 



performed on extremely pure samples having very low compositional disorder (Chandni et al. 20081 have also 



exhibited signatures of stick-slip behavior. While pinning by impurities can never be completely eliminated in real 
systems, one can ask whether such "external" pinning is necessary for stick-slip dynamics in martensites. Can 
there be intrinsically generated disorder, as represented by our <f> field, which induces stick-slip dynamics during 
structural transition in pure systems? 

In order to address this question we need to use our phenomenological theory to study the detailed dynamics 
of the parent -product interface using the droplet ansatz. Calculations in this direction are underway and will be 
published elsewhere. 
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A Appendix: 



The computation of the contribution J^(0o, /3 X , fiy) of the non-afHne field <j) to the free energy of a rectangular 
droplet of size L and W with Af twins in the limit of ai c — > is straight forward, though tedious. The final answer 
is listed here for reference. 



ai<j> 2 + 2aief(f> + c x 



dx 



(d<t>\ 



dxdy 



= [47rai A + 01 Xh c (/3 X , M)p c (J3 y , <j> 

+ ^[aMP^pWvM) + c x h g (p x ,Af)p(p y ,Af) +c y h{p x ) P9 (p y ,N)]<t>l 

where the functions of j3 x , /3 y and M in the square brackets which are all positive definite are, 
h(f3 x ) = 2V^P x {l-e- w2 ^), 



h c (/3 x ,Af) 



h g {l3 x ,N) 



q 2 /3l/2 



2 1- erf(qP x /V2) + e 



-qW 



erf 



V2/3 X 



+ e 



q\V 



erf 



s/2p x 

J [Wl + {W 2 - 2pl) , 



(25) 
(26) 

(27) 
(28) 
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Defining z = q/3 y and u — n/4z, we have further, 

/tt/3. 



p(/3 y ,Af) = \±e- u2 {erf(w) + erf(3u)} + 2{erf(2u) + crf(4w)} + 4 e - 4l,2 erf (2u) 

_ 2 

+ e 2 {erf(2ii + zs) + erf(2ti — iz) — erf(4w + iz) — erf(4tt — ia:)} 

+ 2ie 2 u {erf(3u + iz) — erf(3tt — iz) — erf(u + iz) + erf(u — iz)} , 



(29) 



p c (P y ,N) = [2|erf(\/2M) + crf(2\/2u)| 

+ e" 2z2 | erf (\/2u + i-\/2z) + erf (\/2u - iv^2 

- erf (2V2u + i\/2z) - erf (2y/2u - iV2z) j] , 

P fl (A,,AO = ^ e_4 " 2 " Za [2e z2 {-4(^- 2 2 2 - 4/)erf(2 M ) 

+ e 3 " 2 (-7T 2 + 8z 2 + 16z 4 ) (erf(w) + erf (3m) ) 

+ 4/(1 + 2z 2 )e 4u2 (erf(2u) + erf(4u))} 

+ (8z 2 - 4tt 2 + 42 2 e 4 " 2 ) (erf (2u + iz) + erf(2w - iz)) 

- 4z 2 e 4 " (erf(4u + iz) + crf(4u — iz)) 
+ ie (erf(u + iz) — erf(w — iz)) 

- i(7r 2 -8z 2 )e 3 " 2 (erf(3w + iz)-crf(3u-i2))] (30) 
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